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Abstract 

Invariant and covariant forms of the quark-antiquark interaction derived 
by the method of Eichten and Feinberg are considered. Relations between 
the various terms imposed by Lorentz transformation constraints, here called 
Gromes relations, are found to exist in neither case. Details of the Gromes 
relation proper are reconsidered and inconsistencies found that lead to a vio- 
lation of covariance. 



1 Introduction 



Constituent quark models provide valuable insight into many hadronic phenomena - 
the mass spectrum as well as decay and transition observables. With the appearance 
of advanced probes such as CEBAF, where high resolution data on fine and hyperfine 
spectra to the order of 200 - 300 kev. are expected with complimentary data on 
hadronic production, their predictions are of increasing interest and in increasing 
use. In both sectors non-perturbative QCD effects are expected to be significant. 

Theoretical input for modeling the non-perturbative regime has come from the 
strong coupling limit of lattice gauge theory. Wilson loop area asymptotics |1|] for 
heavy static constituents lead directly to the popular linear confinement model. 
When quarks are given motion their Wilson loop and the potential interaction de- 
rived from it are then manifestly non-local, embodying the confining non-perturbative 
gluon dynamics. A proper accounting for these effects to 0(m^^) has been made in 
the well-known work of Eichten and Feinburg (E.F.), ref|0], where they enter implic- 
itly. They do not however enter into reductions of relativistic equations that input 
only the linearly confining potential; there the semi-relativistic corrections are purely 
kinematical. Hence an increased interest in the E.F. spin dependent interaction and 
the simplifying relation of Gromes 0]. 

It would be difficult to overstate the impact these two results have had on parti- 
cle physics calculations since their appearance and general acceptance over a decade 
ago ( Spires preprint database alone, e.g., lists a combined ~ 700 citation entries.); 
it extends from phenomenological modeling to lattice calculations. The aim of the 
present article is to help clarify their meanings and to set forth criticisms in the 
most accessible terms available. The results are reviewed separately beginning in 
the next section with that of ref[0 where it is demonstrated that whether in its 
invariant or covariant form there are with the potential no accompanying Gromes 
relations(G.R.). In the section that follows the G.R. proper is treated by way of 
its derivation and is found to violate Lorentz covariance. The violation is in con- 
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sequence of two main errors in the derivation: 1) discrepancy between the Lorentz 
transformation and the one employed there as such, and, 2) imposing of Lorentz 
invariance where covariance is required. 



2 E.F. spin-orbit potential 

Beginning from the gauge invariant quark- ant iquark transition amphtude, a hnearly 
confining potential for static quarks is defined in terms of the Wilson loop 

= (0|T*^'=(?/2)P(|/2,2/i)V'(l/i)^(a;i)i'(xi,a;2)V''=(x2)|0) (1) 

^ {tTDtrx Soivu xi; A)P{xu X2)C-^So{x2, 2/2; A)CP{y2, yi)) (2) 

= -tTD ® /) e-^(™^+'"^)^5(xi - yi)5(x2 - y2) , (3) 

6(r) = -^ln(/) = -^\n{trPeM^9jdz,A^{z))) ^ -^ln(l)H^ (4) 

for large T = x° — The symbols above have the following meanings: T* 
time orders, C charge conjugates, the Ps are path-ordered exponentials, P{x, y) = 
Pexp[ig Jy dz^A^{z)], troitrx) is the trace operator for Dirac (gauge) matrices, and 
the average is taken over the gauge fields, {U) = f[dA^]Ue~^'^^''^^\ It is to be 
understood that A^ = A^A" where are the QCD gauge fields and {A'*} are the 
representation matrices for the fermions in the fundamental representation of the 
gauge group SU(3). In (H) Sq is a static fermion propagator. The analysis of ref0 
begins with the introduction of propagators for quarks in motion. They satisfy 

iP-m)Six,y;A) = 5\x - y) (5) 

and may be expanded about the static propagator 

S{x, y- A) = So{x, y;A)+J d'zSo{x, z; A)j ■ BS{z, y- A) (6) 

where Sq obeys 

(L>o7°-^)5o(x,y;A) = 5\x-y) (7) 
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and is given explicitly by 

So{x,y;A) = -^e{x' - y^)e-''-(^'~y'^^^^P{x' ^y^i^ - y) 



(8) 



To order m ^, with > y^, the non-static propagator is given by 



1 + 



ga-B) 



Am? 

+ / d^uS++{x,u;A) 



S++{x,y-A) = S^+{x,y-A) 



2m 4m^ 



(9) 



where the projections are, 5++ = S 



1+7° g 1-7° 



2^2''-' — 2'-'2' ^'^'-^ fo^th. To 

obtain leading relativistic corrections to the static interaction, e(r) of (^, the static 
fermion propagators of (Q) are replaced by these nonstatic ones. Then 

Gi^c = {trntr,{St+ + St- + + Sr) (10) 
X P{y^,y2)C-\St^ + S^- + S,^ + 52-)CP(x2,Xi)) 
, A+7° ^ 1-7% , 1-7° ^ 1-7% , 1 + 7° ^ 1 + 7% ^ 

= -trp{ii^c + Li + ^2) = -trp{ii^c + 4c) 



-fr (7 



(11) 



where Ii^c{Is,c) is identified as the antiquark-charge-conjugated large(small) 2x2 Pauli 
component and tr^ is the trace operator for Pauli matrices. When the large com- 
ponent only is accounted for the combined static and 0(m~^) spin-orbit corrections 
are 



V 



Si ■ Li S2 ■ L2\ e'(r) 



2ml 



r 



+ 



^ 2mf ^nu2 

'Si ■ Li S2 ■ L2\ _^ /S2 ■ Li 
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r \m1m2 mim2 J r 



Sl 



mi T J-T/2 



T/2 



(i2;(z5(B(xi,z))i^/(l),4/ 



2) 



(12) 
(13) 

(14) 

(15) 



where, Lj = r x pj , and 
\ mf 7112 J r V2mf TJ J~t/2 

-(1 ^ 2)) (16) 

V /si 1 /■ rT/2 

(S2-Li-Si-L2)^ = -U--^// rf^rfz'(B(xi,z)D2(x2,^'))W(l)v^ 
r \ 2 I J J-T/2 

-(1 - 2)) . (17) 

A few remarks concerning this result are in order. First, the beginning four 
point function ([^) differs from that of ref where the antifermion fields have not 
been charge conjugated. This difference in four-point functions naturally yields 
a corresponding difference in the interactions derived from them. Second, in the 
interaction potential derivation of [@] there are several algebraic errors leading to 
the final result there (see appendix for details). The interaction resulting from the 
algebraically correct non-charge-conjugated derivation we identify in this paper as 
the E.F. result, Vef of the appendix. This E.F. potential contrasts with the above 
potential most strikingly in that quark and antiquark contributions (for each term 
type) appear in the E.F. potential with identical algebraic sign. In the potential 
above however, derived in the appendix as Vef,c, they appear with opposite sign as 
one would expect for a qq state. For this reason the four-point function of (|I]) is 
seen to have a more acceptable interpretation as a qq propagator than the beginning 
non-charge-conjugated four-point function of ||2| leading to the E.F. result. 

In ref0 it is pointed out that insofar as the derived interaction is related to a 
v.e.v. in a Lorentz invariant theory the interaction itself must also be invariant. 
In fact, V as given in (p!2|), is not Lorentz invariant without a G.R.. It could be 
argued that omission of the small component Is,c in the potential definition (^) 
has ruined the invariance since Is^c appears in the invariant 4-point function (p!T|) 
on an equal footing with the large component. That this is not the case is shown 
by simply replacing Ii^c in (|T^) with Gp^c of (|lTj) and then testing the resulting 
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potential for invariance. The addition to V on this replacement exactly cancels 
line (]TB|), the classical spin-orbit and Thomas precession line, resulting in a spin- 
orbit interaction which we denote Vj^c- Details of the derivation are given in the 
appendix, equations (|77|)-( |101]) . We now consider a Lorentz boost of V/ c to leading 



order in velocity. Care should be taken when transforming the static propagators 
since by construction they must remain solutions to the non-covariant equation (|^). 
To leading order their transformation is 

S^^ix,y;A) ^ S^^{x,y; A) ± J d'zS^^{x, z; A)v ■ z)S^^{z,y; A) (18) 

which is effectively carried out on their time-like path ordered exponentials as 

P(x°,2/°) ^ P(x°,i/°)-z/ dzP{x\z)v ■'D{^,z)P{z,y^) . (19) 

And so the effective transformation (p!9D, carried out in (p^, cancels with the mo- 
mentum transformation of (|1^, and the remaining magnetic field transformation of 
( PI^D yields Vj^c Vi,c + (si ■ [r x v]/mi — S2 ■ [r x v]/m2)e'(r)/r, and noninvariance. 

What must be remembered about the v.e.v. Gj^c of (|T0|) is that its Lorentz 
invariance depends crucially upon its Dirac (or equivalently, its Pauli) trace operator 
without which it is not invariant 

Gi,c ^ Gi,, + 6G (20) 
6G ~ tvp l2,2 ® — ■ / I ci2;(E(x2, z))w X V (21) 

V 1712 "'-f 

+— ■ [ I dz{E{xi, z))wx^® 12x2) ^ (22) 

mi J-f / 

(see appendix for details). From the potential definition (|12|) both as given or with 
Gp c in place of the large component Ii^c it is clear that Lorentz invariance of this 
potential related as it is to Gj^c is not to be expected, cannot be argued as it stands, 
i.e., without the operation of the Pauli trace on J; c ( or Gp^c )• In such an operation 
of course any spin dependence is averaged over, destroyed. Hence there remains 
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no argument for the Lorentz invariance of the spin-dependent contribution to Vj^c 
derived as it is from the invariant Gj^c, and so no argument for Gromes relations 
concerning Vj^c- 

While Vj^c like Vef,c is not Lorentz invariant, it does account consistently for 
the 0(m^^) contributions to the invariant 4-point function from which it is derived, 
while Vef,c does not. Considered as a QQ interaction Hamiltonian however V/ c has 
consistency problems of its own. A simple example will suffice to make the point. It 
is well-known that in the heavy antiquark limit, 1712 — * 00, lines (12) and ([T3|) agree 



exactly with the standard mj"^ expansion of the Dirac equation in Pauli form in the 
presence of a central electric field. Then with the absence of line (|13]) in Vj^c there 
is disagreement between V/ c and the standard reduction of the Dirac equation. 

There is no such disagreement with the interaction derived from a four-point 
function whose Lorentz spinor structure is maintained (see appendix for details). 
This potential we denote Vcov since it is derived from the covariant four-point 
function. Reduction of the covariant four-point function's off diagonal elements 
to 0(m~^) is performed by the Foldy-Wouthysen transformation, U = exp{ts{^)), 
where, s(^) = t'j^'y ■ D(^)/2m with a resulting interaction equivalent to V of 
([1^ . I.e., the covariant interaction is equivalent to the invariant interaction when 
the small component contribution to the invariant trace is (improperly) omitted. 
Its Lorentz transformation properties are unambiguous. From the lowest order 
transformation, pj ^ pj — m^v, B B — v x E, and (|19|) , V transforms as, 
V ^ V + {si ■ [r X v]/2mi — S2 ■ [r x v]/2m2)e'(r)/r = V + V (notice that neither 
Vi nor V2 contribute to ) in strict agreement with its covariant transformation as 
a product of spinor products 

V r-. S++ ^C-^S^'C LS++L-^®LC-^S2'CL-^ (23) 

^ (S++ - iv . 75r+ - i5+-v . 7) ® (24) 

C-\S2- + \S2^^ • 7 + • iSt-)C{2b) 
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~ V + v 



(26) 



where L ^ 1 — a ■ v/2, again, as in the invariant case, precluding the existence of 
Lorentz constraining relations between the transformed terms. 

3 The Gromes method of derivation 

In view of the previous discussions a detailed examination of the G.R. still seems 
worthwhile if only to answer to the enormous interest it has generated over the last 
decade and a half. It should be clear however that such an examination is little 
more than an academic exercise. In refQ] the lowest order Lorentz transformation 
on V by which the relation is derived is implemented in three distinct steps. 

1) momentum transformation: 

Pi Pi - rriiW . 

2) magnetic field transformation: 

B ^ B - V X E 

which by the given lemma, (2.8) of @, is effectively carried out on the relative 
coordinate as 

r ^ r + V X Si/mi — v x S2/m2 . 

3) Tj transformation: 

^ Tj + [v X Sj]/2mi 

following the definition and transformation specification of the quantity 

X = Ti + {pi/mi)t + [pi X Sj]/2mi 
X X — vt . 
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Effects of transformations 1) and 2) taken together with the effective transforma- 
tion of the path-ordered exponentials, (JT^), combine in V of (0) to give the correct 
covariant transformation (|23|) , V —>■ V + V. 

Transformation 3) however is entirely spurious with no correspondence to the 
field theoretic Lorentz transformation. Its effect is to produce the extra transfor- 
mation terms, e(r) —>■ e(r) + ([v x Si] ■ r/2mi — [v x S2] ■ r / 21712) e'{r)/r. The G.R. 
appears when 1), 2), and 3) are taken together (with a sign error included for the 
terms from step 2)) and invariance imposed to 0(m^^). Then, e + Vi — V2 = 0. 

The relation has received derivation in several other contexts as well. These 
derivations should each one undergo equal scrutiny. 



4 Summary 

The E.F. interaction Hamiltonian contains dynamical information on the nonper- 
turbative gluon field not present in models that simply relativize the static linear 
potential. But in its given form one is left either to make simplifying assumptions 
or to pursue evaluation on the lattice. For the latter, it is unclear whether the sign 
discrepancies between the results here, ([T^) - ([T6|), and ( Vi, V2 ) of ref0 would 



affect corresponding discrepancies between the respective lattice reductions. For 
the former, the electric confinement ansatz, Vi = V2 = 0, in conceptual agreement 
with Buchmiiller's picture^], assumed in and later abandoned (due to variance 
with spin phenomenology), demonstrates that these should be made only with spe- 
cial care and that they may in any event have consequences difficult to predict. In 
these regards it may be viewed as a shortcoming of the E.F. formalism that here too 
only the static limit of the minimal area law enters explicitly; nonlocality is implied 
only (e.g., via electric field insertions 0). This happens as an artifact of the initial 
fermion propagator expansion (^ around the static limit, and is readily remedied 
when these propagators are replaced by ones more compatible with semi-relativistic 
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fermion motion. Then the Wilson loop of (^) is explicitly nonlocal and relates the 
resulting potential directly to the minimum area. 

Such a program has been carried through in an article by Brambilla, Consoli, 
and Prosperi |^ both for spin and spin- independent corrections to 0{m~^). From 
propagators expressed as integrals over the phase space and the Nambu-Goto ac- 
tion as the effective area their final potential is given in terms of familiar quantum 
mechanical operators (It should be pointed out that the spin-orbit agreement with 
ref0 follows from a systematic mathematical error made in their appendix ||^. See 
also ref0. ). 

A more rigorous implementation of the electric confinement ansatz might also 
be attempted, though it is uncertain whether Buchmiiller's picture remains self- 
consistent for quarks in motion or whether it agrees with minimal area asymptotics. 
These questions are considered in reference 

Acknowledgments: This work was supported in part by the National Science Foun- 
dation under Grant No. HRD-9154080. I am also deeply indebted to The College 
of William & Mary and in particular to F. Gross for providing me office space and 
full access to the Physics department facilities during the course of this work. 

5 appendix: potential derivations 

In evaluating the gauge-invariant QQ 4-point function 

G = {0\T*^{y2)P{y2,yi)^{yMxi)P{xi,X2)^{x2m (27) 
{trS{x2,y2;A)P{y2,yi)S{yi,xi;A)P{xi,X2)) (28) 

the trace appearing in (p8[) may be taken over both gauge and Dirac matrices or 
over gauge matrices only, depending, respectively, on whether the four suppressed 
Dirac indices of (^) are of two like pairs or are all distinct. In the first instance the 
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four point function is invariant under Lorentz transformations and in the second it 
is covariant, transforming as a product of spinor products. Interaction potentials 
derived from these invariant and covariant 4-point functions therefore have trans- 
formation properties that are invariant and covariant, respectively. Both 4-point 
functions appear in the literature as starting points for the potential derivation; for 
example, in the invariant version is used, and in the covariant version. 

Interaction potentials from the invariant 4-point function : 

First we derive interaction potentials from the invariant 4-point function. The 
derivations follow closely that given in pf . Fermion propagators of (pSD are expanded 
around the static propagator as 

S{x, y; A) = So{x, y; A) + J d^zSo{x, z; A)-f ■ BSiz, y; A) (29) 

with 

So{x,y;A) = -^e{x' - yy-^"^^^'-y'^^^P{x' ,y')6{^ - y) 

-le^y^ - xO)e-*'"(^"-"")^^P(x°, y°)5(x - y) . (30) 

They are then expanded over static energy projectors 

S = l + 7%l+7° l + 7° ^ l-7° l-7° g l+7° l-7° g l-7° 
22222222 

= S++ + S+- + S-^ + S— (31) 

so that 

5++(x,i/;A) = S^+ix,y;A) + J d'zS^^ix,z;A)j-Ti{z)S-'-iz,y;A) (32) 

S-'-{x,y;A) = J d^zS++ (x, z; A)-f ■ Biz)S- iz,y; A) (33) 

S-+ix,y;A) = J d'zS,-ix,z;A)-fBiz)S++iz,y;A) (34) 

S~-{x,y;A) = S,- {x,y; A) + f d'^zS,' {x, z; A)-f -D{z)S+- {z,y- A) (35) 
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The projected propagator components above are evaluated to 0(m~^) by substitu- 
tion of the static propagator solution ( ppp and iteration to this order for x° > y^. 
For the 5*"^^ component then 

5++(a;,t/) = S^+{x,y) + J d^wA-^~{x,w)j ■'D{w)S^+{w,y) (36) 

where 

1 + 7° 



-e 



2 

1+7° -im{x°+w'' 



-e 



)5(x-w)/ (37) 



and 

/(z°) = P(x°,z°)7-D(z)P(2°,w;°). (38) 
We rewrite the above integral as 

/ = e2''""'"^(x° - u;°) r dzV'"^V(^° + w°) (39) 

J— oo 
v/ — oo 

and use the 0(m~^) integral relation 

/;^..v-7(.»+o - j^/'tt) (41) 

with 

/(ti;°) = P(a;°,^°)7-D(w;), /'(«;°) = -zP(x°, ^°)[Z}°(u;), 7 ■ D(«;)] (42) 

/(x°) = 7-D(x)P(x°,w;°), /'(x°) = -«[D°(x),7-D(x)]P(x°,t/7°) (43) 

to find on substitution 

A+-(x,^) = -S^+{xM{^^ ■Ti{w) + ^[D\w)n --Diw)]) (44) 

-(^7 ■ D(x) + ±-^[D\x)n ■ D(a:)])5o-(x, 1/;) (45) 
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and 



1 



■[D%w),-f ■ DH])7 • Biw)S++{w, y) 



where 



Aw? 

-(^7 ■ D(a;) + -^[D\x)n ■ D(x)]) A-+(x, y) 



A-+{x,y) = J d^wSo'{x,w)-f'D{w)S^+{w,y) 



2m 



With this, and from the identities 



(7-D) 



j--D{x)S++{x,y). 



+ g(TB 



D°,7-D 7-D 



-ig{6ij - ^e^jkcr^)E'D^ 



we arrive at 
1 



1 + 



Am? 



S++{x,y) = S^+{x,y) 



(D' -^cr-B) 



+ y (i ujSq {x, ijj 
equation (H). Then finally, to 0(m^^) 



^9 
Avn? 



— (D^-oCT-B) 



Am? 



Am^ 
1 



Am^ 



^ / d^ud^zS^+ix,uj)iTi^ -g(T-B) 
J 

xS++{uj,z){B'-g(T-B)S++{z,y) 
B'-g(T-B)S^+ix,y). 
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By similar methods we find ( x° > y° ) 



and 



S+-{x,y) 
S-+{x,y) 
S-{x,y) 



S ^{y,x) 
S+-{y,x) 
S^+{y,x) 



-^S^^{x,y)-i-Ti{y) 

-^7-D(x)5o++(a;,|/) 



1 



Am? 



Am? 



7-D(a;)5o++(a;,y)7-D(y) 



~S^-{x,y)[J^\y)-ga-B 



+ig I d''zSo-{x,z){6ij-i€ijkcr'')E'{z)So-{z,y)D'{y) 



S {y,x) ~ {y,x) + j (I'^ujSq {y,cj) 



2m 



^9 



Am? 



Sq {uj, x) 



Arn? 
1 



^ j d'ud^zS^-{y,u){^^ -g(T-B) 

xS^-{u,z){p^ - ga ■B)S^-{z,x) 
^^-ga-B)S^-{y,x) 



Am? 

-7;-SQ-{y,x)-i ■T>{x) 
2m 

-2^7-D(y)5o"(y,^) 
1 



4m^ 



-r-T^{y)S,~{y,x)-fT>{x) 



Am^ 



St^{y,x)[T^\x)-ga-B{x)] 

-ig / d''zS^^{y,z){5ii-ieiika'')E\z)S++{z,x)D^{x 

where the identity 

-i-'D{x)S^^{x,y) = 4^^(x,t/)7-D(t/)±z|ciS4^^(x,^)7-E(^)5, 
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with 



oa-7° 



P(x",i/")5(x-y) 



has been used, and the off- diagonals, and S + are shown to lowest order. Note 
that in evaluating S {y-,x) above the integral relation used in place of equation 



(^iD ( which is used in the evaluation of S^^{x^y) ) is 



-/(O 



(76) 



/o 2m'' Am?' 

yielding the algebraic sign difference between lines ( [55| ) and (]65|). For construction 
of the potential interaction we express the invariant Green's function as 



Gi = {trSiPS2P) 

= trnitr^S.PS^P) 

= -trBGe-*('"^+™^)(^"-J^°)(5(xi - yi)5(x2 - y2) 

= -trpGpe-^(™^+™^)(^°-^")5(xi - yi)5(x2 - y^) 



(77) 
(78) 
(79) 
(80) 



where in ( [78D Dirac and gauge traces have been explicitly separated, tvp is the trace 
operator for Pauli matrices, and Gp is the sum of 2x2 Pauli components along the 
diagonal of G. The "invariant" potential interaction is then given by 



Vr 



(81) 



for large T. Substitution of 
function 



into (2S) then yields for the invariant 4-point 



Gj = {trDtrx[S^^{xi,yi) + S (xi, ?/i)]P(t/i, t/s) 

X (^2, X2) + 5++(y2, X2)]Pix2, Xi)) 

= -troitrxP 



4-^0 1-^0 ■ 

2 "-^ip ' 2 '^P 



2 ^2p 2 2p 



(82) 

(83) 
(84) 



X exp(z(?^d^^A'^(^)))e-'('"i+"^^)^5(xi - yi)5(x2 - y2) (85) 
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where x° = —y^ = ^ has been taken, and the "p" subscript identifies Pauh compo- 
nents which are given by 



9+ 



1 - 



2mi 



dz(D'^{xi, z) — (Ti ■ B(xi, z)) 



Am{ J-'L 



dz eijk(y\E\^i,z)D^{^i,z) 



S- 



S. 



2p 



4m\ 
1 



dz eijkCTfE\^,,z)D\^,,-T/2) 



2m2 



(i2:(D^(x2, z) - (T2- B(x2, z)) 



^2p 



Am? 
i 



2 

2 

T 



dz eijkcr^E' (xa , 2;) (xa , 2;) 



Ami 



dze,,kalE\^2,z)D\^2,T/2) 



(86) 
(87) 
(88) 
(89) 
(90) 
(91) 



2 •'-T 



for 0(m~^) static and spin-orbit contributions only. Gauge fields evaluated at 
have been set to zero. 

To show explicitly that the Lorentz invariance of Gj depends upon the operation 
of its Pauli trace we consider the leading order boost L^l — Q;-v/2onG/ 



1 , 



Gi ^ Gj + SGj = {tr Si P S2 P) 

® P ([5-+, V • 7] - [5+-, V ■j])P + 

P([5-+,v7]-[5+-,v7])P®5o- 
giving on substitution for S*^^ and the use of identity (ff^ ) 

6G 



tvp l2x2 ® I dz{E{x2, z))w X V 

^ m2 



/ ^ dz{E{yii, z))w X V (g) 12^,2 

mi 



0, 



(92) 
(93) 

) (94) 

(95) 
(96) 



From fBOD then 



Gr, 



{{St + Si,) ® (^2-. + st)w) 



2p "T ^2py' 



(97) 
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^ h + Is (98) 

with 

= {St,®S^^W) (99) 

/. = {iS^, + Sl)W) (100) 

where J; is the large component ( identified by its nonzero value in the static limit, 
m — > oo ) and W = exp{zg § dz^A^{z)) is the Wilson loop. From ( [HT| ) we have for 
the interaction 

Vi = -^\n{Ii + Is) (101) 
e(r) (102) 

/ Si ■ Li 82" 1-12 \ V( / S2 ■ Li Si • L2 \ V^' 



mf / r \m1m2 mim2 J r 



(103) 



^ Si 1 /•^/2 \ 

— ■- dz{zgB{^„ z))w/{l)w + (1 - 2) (104) 
^mi I J-T/2 J 



where, Lj = r x pj , and 
/si -Li S2 ■L2\ V; _ / si If rT/2 



j ^ dzdz'{B{^^, ^)D2(xi, z'))w/{l)w 



m\ 7712 J \2mi TJ J-t/2 

+ (1^2)) (105) 

V' /Si If /•^/2 



V /Si 1 /■ ^ 

(S2-Li-Si-L2)^ ^ -(i'Ti y^//^^^'(B(xi,^)D2(x2,z'))W(l)v^ 

+ (1 ^ 2)) (106) 

and the relation 

I I' dzei^k<y\E\^,z)D^{^,i)) -Te-^M^cr-L^ (107) 



has been used. The sum // + Is appearing above in (|9q) is to be compared with I 
in equation (4.3) of p[. When the "small" component of the invariant trace, J^, is 
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ignored Vj becomes 



Vef = e(r) (108) 

y mf m2 ) r \m1m2 mim2 J r 

- ^ ■ - / rf^(2^?B(xi, z))w/{l)w + (1 2) (111) 

with 

f^-^)^ = -(-^.l/r^.rf.'(B(x„.)D^(x„/)W(lV 
y mf 7112 J r \2mf TJ J-r/2 

+ (1^2)) (112) 

T/' /Si 1 r /■2"/2 

(S2-Li-Si-L2)^ = -hr-^// rf^rf^'(B(xi,z)D2(x2,^'))W(l)iy 
r \ 2 1 J J-T/2 

+ (1^2)^. (113) 

The "large component only" interaction is precisely what is derived by authors 
Eichten and Feinberg in 0. It's differences with the above Vef potential are easily 
traced to three algebraic errors made in their derivation: 1) In their appendix there 
is an expression for S^^, (A3), given in terms of an integral whose evaluation is 
given in (A15). Direct substitution of (A15) into (A3) does not yield their equation 
(2.11) for S^^, but yields instead equation ( [53| ) above with spin terms on the rhs 
differing in algebraic sign. 2) There is a misrepresentation of the coordinate space 
momentum operator in going from equation (4.11a) to (4.11b) in the E.F. work. 3) 
There is a missed sign in the evaluation of an integral for S {y, x) whose analogue 
in the evaluation of S^~^{x,y) is (A7b). Here this mistake amounts to using the 
integral result of (^1]) in the evaluation of S {y, x) where the result of ( [76D is called 
for. These individual errors would have the following effects on the above Vef'- i) 
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error 1) would change the overall algebraic signs of lines ( |109| ) and( |lll| ). ii) error 2) 



would partially correct error 1) by again changing the algebraic sign of line ( [I09| ). 
iii) error 3) would change the relative sign between terms in line ( |109| ) making the 
subscripted "2" term negative. The resulting interaction Hamiltonian would then 
be identical with the static plus spin-orbit interaction of reference 0. 

The conspicuous lack of relative sign differences between subscripted "1" and "2" 
contributions in Vef above makes its interpretation as a QQ interaction Hamiltonian 
problematic. When the antiquark field operators in the beginning four point function 
(pTj) are properly charge conjugated ( as they are for example in |0)this problem 
disappears. Specifically, the necessary change is 

^{X2) ^ r{x2) (114) 



T 



^S{x2,y2;A) C-'S{x2,y2;A)C = - [S{y2,X2;-A')\ (115) 
where the outer transposition in the last step is taken over both gauge and Dirac 



matrices. The effect on Vef of ( |108| ) from the above charge conjugation made in 
(p8D is simply to change the algebraic signs of the field insertions on antiquark lines 
of the Wilson loop while preserving the original path ordering 



2m? 2m2 J r 



Vef,c = e(r) (116) 
'si ■ Li S2 ■ e'(r 

y mf 1712 / ^ \m1m2 mim2 / r 

/ Si 1 /-^/^ \ 
- — ■ - / dz{tgBi^^, z))w/{l)w -(1^2) (119) 

with 

/si -Li S2 ■L2\ VI _ / si If 



mf 7712 I \2m 



1 r ni'^ 

2 ■ - / / dzdz\B{^^, ^)D2(xi, 2!y)wl{\)w 

-(1 ^ 2)) (120) 
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V /Si 1 /■ r^/^ 

(S2-Li-Si-L2)^ = -(^-j--J J ^^^dzdz'{Bi^^,z)B'i^2,z'))w/{l)w 

-(1^2)). (121) 

where the "c" subscript indicates that antiquark fields have been charge conjugated. 
This is our V of equation (p^. 

To obtain the full Vj charge conjugated interaction likewise simply requires that 
we change the algebraic signs of the field insertions along the antiquark lines ap- 
pearing in Vj above 

Vi,c = e(r) (122) 



+ 



nil ^2 j \m1m2 17111712 J r 

■ ^ dz{tgB{^i, z))w/{l)w - (1 - 2)) (124) 



with 

/Si-Li S2-L2\ VI _ / Si If 



1 r ri"^ 

I'tI j_^^/''dz'{B{^i,z)J^\^i,z'))w/{l)w 
-(1 ^ 2)) (125) 



S2-Li-Si-L2)^ = -{i-7f dzdz'{B{^,,z)-D\^2,z'))w/{l)w 
r \ 2 I J J-T/2 

-(1 ^ 2)^ . (126) 



Interaction potential from the covariant 4-point function : 

Here we derive the interaction potential from the covariant 4-point function 
whose antiquark fields are charge conjugated. 

Gcov = {trxSixi,yuA)P{yi,y2)C-'Siy2,X2;A)CPix2,Xi)) (127) 
_ _G'e-(™i+-2)^5(xi - yi)5(x2 - ya) . (128) 
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The interaction is found from the large component of the diagonahzed G by 

Ko.,c = -^InG, (129) 

where Gi is the large component. The diagonalization to 0(m^^) is performed via 
the Foldy-Wouthysen transformation U = exp(zs({)), with, = I'f^'y ■ D({)/2m. 
It will suffice to reduce the off diagonal elements to 0(m~^), accomplished by a 
single transformation. For the quark propagator reduction, using the above results 

S++(x,y) = S++(x,y) + ld'wS+-^(x,w)^{^-T>{w))'S++(w,y) (130) 
+0(m-2) (131) 
S^-ix,y) = -l-S^+{x,y)j-Biy) + Oim-') (132) 

S-+{x,y) = -l-j.B(x)S^+(x,y)+0(m-') (133) 

S—{x,y) ~ 0(m-2) (134) 

we have 

S{x,y) 5'(x,y) = e'"(")5(a;,i/)e-^"(^) (135) 

= S{x, y) + i[u{x) (^++(x, y) + 5+-(x, y) + 5-+(x, y)) (136) 

-(5++(x, y) + 5+-(x, y) + S-+{x, y))u{y) (137) 

-^Kx)(i.(x)5++(x,y) - S++{x,y)u{y)) (138) 

-(m(x)^++(x, y) - ^++(x, |/)m(i/))m(|/)] + 0(m-3) (139) 
= '^^^(^'^) + ^[(^ • ^{x)fsr{x,y) + ^o^+(x,y)(7 • ^{y)f] 

+ ^ J d'w[S+-^(x,w)(-fB(w)yS++(w,y)j-T>{y) (140) 

-7 ■ D(x)5(j"+(x, w)(7 ■ T){w)yS^+{w, y)] + 0{m-^) (141) 

The off diagonals are now of 0(m~^). Another transformation would further reduce 
the off diagonals to 0{m~^) while contributing nothing to the diagonal. The diag- 
onahzed quark large component is therefore given by S'^'^{x,y) alone ( with gauge 
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fields again set to zero at ±^ ). The antiquark propagator diagonalization follows 
along the same hnes, leading to the large antiquark component C^^S {y,x)C. 
Then from (jS^) to (|99|) above we find 

Gi = ii,c (142) 

where the "c" subscript indicates that field insertions on the antiquark line of the 
Wilson loop in J/ of (^) have undergone an algebraic sign change. This gives 

Vcov,c = Vef,c (143) 
= V (144) 

where Vef,c is given in ( p.l6|) and V in (|12|). I.e. the interaction obtained from the 
diagonalized covariant 4-point function is identical with the one obtained from the 
large component of the invariant interaction. 
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